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Abstract 

Understanding tlie roles of disorder and metal/grapliene interface on tlie electronic and transport properties of 
graphene-based systems is crucial for a consistent analysis of the data deriving from experimental measurements. 
The present work is devoted to the detailed study of graphene nanoribbon systems by means of self-consistent 
quantum transport calculations. The computational formalism is based on a coupled Schrodinger/Poisson approach 
that respects both chemistry and electrostatics, applied to pure/defected graphene nanoribbons (ideally or end- 
contacted by various fee metals). We theoretically characterize the formation of metal-graphene junctions as well 
as the effects of backscattering due to the presence of vacancies and impurities. Our results evidence that disorder 
can infer significant alterations on the conduction process, giving rise to mobility gaps in the conductance 
distribution. Moreover, we show the importance of metal-graphene coupling that gives rise to doping-related 
phenomena and a degradation of conductance quantization characteristics. 



Introduction 

Graphene nanoribbons (GNRs) are the most promising gra- 
phene-based nanostructures for electronic appUcations 
since they are potentially suited for band-gap engineering, 
maintaining the excellent electronic properties of the parent 
two-dimensional graphene layer. GNRs have been already 
synthesized by means of different pattering techniques [1,2], 
and there exists convincing evidence that their electronic 
structure manifests subband formation which is a typical 
predicted signature of the one-dimensional (ID) confine- 
ment [3]. Actually, a useful intrinsic (i.e., due to the geo- 
metric confinement) band-gap value can be obtained only 
for GNR structures with widths of approximately 1 nm that 
cannot be easily fabricated with the current lithographic 
techniques. As a consequence, the controlled inclusion of 
defects and/or impurities in GNRs has been proposed to 
overcome the minimum of conductivity problem [4]. Any- 
how, we should assume that unintentionally defected GNRs 
are present from the phase of production due to the impos- 
sibility of an accurate control of purity at an atomic level. 

Electrical characterization studies of GNRs must preli- 
minarily consider the impact of disorder in their 
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electrical properties, which is crucial also due the possi- 
ble Anderson localization phenomena in these quasi- ID 
systems [5]. Moreover, the formation of a junction 
between a relatively large metal probes (i.e., a three- 
dimensional (3D) system) and GNRs (a ID system) is at 
the basis of any electrical measurements. Understanding 
the role of interface bonding and electrostatics in the 
contact region is also crucial in order to categorize the 
electronic characteristics of these systems. Indeed, we 
could argue that in the absence of defects, a significant 
source of resistivity in GNRs and the consequent devia- 
tions from the ideal behavior should be derived from 
the interaction with the metallic electrodes. In this 
work, these issues are theoretically investigated and 
their implications on the interpretation of electrical 
measurements on GNR-based systems are discussed. 

Theoretical approach 

We consider pure and defected hydrogen-terminated arm- 
chair GNRs (AGNR) of different widths and lengths (from 
few nanometers to approximately a micrometer) with two 
different contact configurations: (a) ideally contacted 
(same width as the conductor without defect or impurity 
inclusions) at the right side and the left side, and (b) end- 
contacted by 3D semi-infinite metallic electrodes (Au, Pd, 
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and Al (111) surfaces) at the left side whereas ideally con- 
tacted at the right side. The terminology of Ref. [6] is 
applied to categorize them on the basis of the dimer lines 
A/^a along the ribbon width (e.g., with = 45 AGNR, we 
indicate a ribbon with 45 dimer lines). We use the non- 
equilibrium Green function (NEGF) 



G = [E^I - H - El - Ei^) 



(1) 



and the standard Landauer-Buttiker approach for the 
calculation of quantum transport [7]. Hamiltonians Hq 
and overlap matrices r are written within first-principles 
parameterized models [8]. 

Metal surface Green functions for the evaluation of 
the respective self-energies E = x g^T^ are calculated for 
the 3D semi-infinite contact with a back and forth real 
to /c-space Fourier transform exploiting lattice periodi- 
city [8]. Mobile charges Pf deriving from the NEGF are 
passed to a 3D numerical Poisson solver from where the 
self-consistent potentials are calculated by solving V^U^c 
= -pf/s, where H = Hq + Use Equation 1. Electrical 
potential has been fixed at the metal-GNR interface 
(applying a Dirichelet type boundary condition) at a 
value Z/ieft sc = 0m - 0gr> where 0^ and (p^^ are the 
experimentally measured work functions for (111) 
metallic surfaces and graphene [8], while zero electric 
field (i.e., null Neumann boundary condition) has been 
set for the ideal contacts. Self-consistency is enhanced 
by a predictor/corrector Newton-Rapson algorithm. 
Fermi-Dirac statistics have been introduced in the simu- 
lation scheme for temperature T = 300 K. 

Defected GNRs 

Isolated defect and impurity alter both the density of 
states distribution DOS(£) than conductance g{E) of 
quasi- ID GNRs. In particular n-type (p-type) impurities 
introduce electronic states (also called resonance states 
[3]) at energies above (below) the charge neutrality 
point (i.e., the Dirac point) of the pristine pure systems. 
In correspondence of these states, a decrease of the con- 
ductance with respect the pure system can occur 
because of the scattering of the conducting electron by 
the effective potential due to the impurity presence [3]. 
The conductance alteration strongly depends on the 
position of the scattering center in the nanoribbon 
structure [8]. However, in general, we expect that defec- 
tive GNRs are characterized by larger density of the 
states and smaller conductance with respect to the ide- 
ally pure GNRs. 

In a real defected system, we expect a finite density of 
random distributed scattering centers and, as a conse- 
quence, the effect of multiple scattering processes 
should be evaluated in this kind of configuration. In 
Figure 1, we show the (small bias) average conductance 
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Figure 1 Average conductance for a nitrogen-doped = 45 
AGNR. Average conductance <g> as a function of the energy E for 
a nitrogen-doped A/a = 45 AGNR of different lengtlis. Plotted values 
represent statistical averages over more than 500 equivalent replicas 
of the system. Charge neutrality points of pure and doped systems 
are aligned at ^ = 0 eV in the figure. 



of an ideally contacted A^^ = 45 AGNR doped with a 
0.2% density of nitrogen atoms and for systems with 
increasing length L, from approximately 0.1 [im to 
approximately 0.8 [im. An asymmetric decrease of the 
average conductance due to the impurity scattering can 
be observed for the whole spectrum. This behavior is 
particularly important in the energy region near the sin- 
gle impurity resonance states (i.e., for energies E ~ 0.2 
eV) where a mobility gap appears also for the smaller 
systems. We note that the pure N^, = 45 AGNR is a 
semiconductor GNR with an energy gap of approxi- 
mately 0.2 eV, 

A qualitatively similar behavior is shown by the 
vacancy-damaged and nitrogen-doped semimetallic = 
47 AGNR (see Figure 2). In the vacancy-damaged case, 
a large mobility gap appears in the negative energies 
(hole band) region [9] due to the strong backscattering 
of the defects. However, apart from the intensity of the 
scattering, vacancy-defected systems have a p impurity- 
like behavior. 

From an analysis of the conductance spectra, we can 
derive same general features which can be useful for the 
interpretation of the electrical characterization of real 
GNRs. Indeed, we note that the conduction spectrum 
can be measured in a three terminal configuration [4] 
tuning the gate potential in order to modify the electron 
density in the nanostructure. Firstly, we note the persis- 
tence of the conduction modulation with energy in dis- 
ordered systems, which is a marker of the conductance 
plateaus of the subband structure in pure GNRs. The 
effective transmission in the subbands is strongly 
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Figure 2 Average conductance for nitrogen-doped 
vacancy-damaged AGNRs. Average conductance <g> as a 
function of the energy E for a nitrogen-doped = 45 AGNR (solid 
line), nitrogen-doped A/a = 47 AGNR (dashed line), and vacancy- 
damaged A/a = 47 AGNR (point) with fixed length: L ~ 0.21 [\m. 
Plotted values represent statistical averages over more of 500 
equivalent replicas of the system. Charge neutrality points of pure 
and defected systems are aligned at ^ = 0 eV in the figure. 
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Figure 3 Conductance fluctuation for nitrogen-doped and 
vacancy-damaged AGNRs. Conductance fluctuation cr as a 
function of the energy E for a nitrogen-doped A/a = 45 AGNR (solid 
line), a nitrogen-doped = 47 AGNR (dashed line), and a vacancy- 
damaged A/a = 47 AGNR (point) with fixed length: /.~0.21 |jm. 
Plotted values represent statistical averages over more of 500 
equivalent replicas of the system. Charge neutrality points of pure 
and defected systems are aligned at ^ = 0 in the figure. 



reduced due to localization effects that suppress the 
conductance transparency in the spectral region where 
resonance states are located. Finally, the general occur- 
rence of mobility gaps (see also Ref. [10]) can be hardly 
distinguished from the intrinsic bandgaps in semicon- 
ductor GNRs when experimentally measuring the con- 
ductance. A clear signature of the mobility gap is the 
occurrence of huge values of conductance fluctuations 
in the same spectral region of the gap due to its back- 
scattering origin. In order to demonstrate this assump- 
tion, we have plotted in Figure 3 the statistically 
evaluated conductance variance 

=< (In^)^ > — < ln^>^ 

(where g is in 2e^/h units) for the same cases of Figure 2. 
In the quasi-metal systems, the conduction depletions are 
systematically related to variance peaks, while the variance 
in the standard gap of the semiconductor GNRs does not 
show particular features. 

Metal-GNR junction 

Self-consistent local density of the state (LDOS) calcula- 
tions have shown a non-conventional scenario for the char- 
ging and the related electrostatics of the heterojunction 
between the 3D metal and the ID GNRs [8]. In particular: 

(a) The electrostatic potential, due to the work func- 
tion difference, drops in a few-nanometer region 



near the interface and has non-zero flat value, 
denoting the presence of carrier accumulation 
throughout the GNR length; 

(b) Unlike the standard metal semiconductor inter- 
face, band bending is not rigid for both conduction 
and valence bands as a result of quantum interfer- 
ence (between wave functions with ID character in 
the nanoribbon and 3D character in the metal) and 
electron confinement. For example, in the case of a 
metal with high (low) work function with respect to 
the graphene, like Au (Al), conduction (valence) 
band states shift smoothly following the electrostatic 
potential, while discrete localized states appear with 
a few-nanometer spatial extension in the region 
where the valence (conduction) band maximally 
bends; 

(c) Quantum interference, localized states, and the 
metal induced gap states (i.e., the tails of the metal 
wavefunctions) characterize the LDOS in the junc- 
tion region. 

The local electronic structure characteristics in the 
contact region can non-trivially influence the conduc- 
tion mechanism since, e.g., localized states do not con- 
tribute to the conduction, giving rise to conductance 
asymmetries and an overall loss of the transport infor- 
mation with respect to the ideal case. In Figure 4, small- 
and finite-bias conductance spectra of a pure = 16 
AGNR end-contacted with Au are plotted as a function 
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Figure 4 Small- and finite-bias conductance spectra. Small-bias 
(points) and finite-bias (+0.5 V green line, -0.5 red line) 
conductances of a pure /Vg = 16 AGNR end-contacted with Au. 



of energy. These spectra show similar characteristics 
with the ones obtained considering metals with higher 
work functions than graphene (e.g., Pt, Pd). A Schottky 
barrier of the order of 0.2 to 0.3 eV can be determined 
by the difference between the gap in the conductance 
spectrum for the contacted GNRs at 0 V (Figure 4 
dotted line) and the gap (approximately 0.6 eV) of the 
ideal non-contacted GNR. 

Apart from the Schottky barriers, we note that conduc- 
tion band charge flow is strongly suppressed due to the 
band bending and the p-type character, giving an asym- 
metric form to the overall conductance distribution 
(asymmetric spectra, not shown here, are calculated also 
for quasi-metallic GNRs ). Note that finite-bias conduc- 
tance spectra are affected by the presence of the Schottky 
barrier. Indeed, while the -0.5- V case shows the alteration 
with respect to the small-bias case due to non-equili- 
brium charging (see Ref. [11] for a complete discussion), 
the +0.5 is almost identical to the small-bias case. 

Current-voltage (/-V) characteristics of the junction in 
the case of a pure = 16 AGNR contacted with differ- 
ent metals are reported in Figure 5. Larger / values 
obtained for a negative bias in the case of Pd with respect 
Au are due to its slighter more pronounced p-type char- 
acter (while in turn Au seems transparent near the Fermi 
level with the conductance arriving at the 1 Gq = 2e^lh 
plateau of the ideal case). Al has a lower work function 
with respect to graphene and the Al-GNR junction 
shows a quasi-ambipolar Schottky behavior (i.e., the I-V 
characteristic is almost symmetric for positive and nega- 
tive bias). However, in the latter case, the dominant 
aspect is the strong scattering by the contacts and the 
related suppression of the contact transparency. 
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Figure 5 Current-voltage characteristics. l-V characteristics 
derived by the NEGF in the Landauer-Buttiker scheme for a A/a = 16 
AGNR, end-contacted with three different metals Au (black line), Pd 
(green line), and Al (red line). 



Conclusion 

The results of this theoretical study have shown that 
defects, doping, and metal/graphene interfaces can non- 
trivially influence the conduction mechanism of gra- 
phene nanostructures, giving rise to mode transmission 
suppression, conductance asymmetries, mobility gap for- 
mations, and strong conductance fluctuations. The 
implication of these results for the electrical characteri- 
zation of real structures has been discussed. A general 
crucial issue, which should also be taken into account 
when probing transport in GNRs, is that the electro- 
statics and the chemical bonding aspects can act com- 
plementary for the determination of transparency in 
graphene-based nanostructures. 
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